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We compute the spectral function of a spin S hole injected into a two-dimensional antiferromagnet 
or superconductor in the vicinity of a magnetic quantum critical point. We show that, near van 
Hove singularities, the problem maps onto that of a static vacancy carrying excess spin S. The 
hole creation operator is characterized by a new boundary anomalous dimension and a vanishing 
quasiparticle residue at the critical point. We discuss possible relevance to photoemission spectra 
of cuprate superconductors near the anti-nodal points. 



A vigorous debate has been stimulated by photoe- 
mission experiments 0] showing that the gapped anti- 
nodal quasiparticle excitations of the high temperature 
superconductors have a broad energy distribution curve 
(EDC). The many proposals for this anomalous behavior 
include (i) electron scattering by antiferromagnetic fluc- 
tuations ||, (ii) one-dimensional (Id) fluctuations at in- 
termediate scales leading to Id electron fractionalization 
(Hi) 2d electron fractionalization induced by proxim- 
ity to exotic 2d spin liquid states p|, and (iv) coupling 
to superconducting phase and vortex fluctuations || . 

The studies of this paper fall into category (i); how- 
ever, our results also extend into a quasi Id regime, and 
offer a different perspective on category (ii). We will 
examine the spectral function (or EDC) of the gapped 
anti-nodal quasiparticles [|| in the vicinity of a quan- 
tum critical point between a d-wave superconductor and 
a state with co-existing superconducting and antiferro- 
magnetic order. Our results will also apply to insulators 
in the vicinity of a quantum transition between states 
with antiferromagnetic order and a spin gap: the spec- 
tral functions are then those of holes (and electrons) near 
van Hove singularities in the band structure of the spin 
gap state. Our perspective on the superconducting case 
differs from that of [|| in that we depart from ground 
states in which pairing correlations have already induced 
a gap to fermionic excitations in the anti-nodal regions, 
rather than from states with a Fermi surface. Appealing 
to the proximity of a magnetic quantum critical point 
between the former states allows us to make controlled 
statements in a regime with strong coupling between the 
fermionic quasiparticles (henceforth referred to generi- 
cally as holes) and the antiferromagnetic fluctuations. 

Our primary result is that the proper characterization 
of the damping of the holes by low energy spin fluctua- 
tions is provided by the framework of "boundary" critical 
phenomena |7j]. The hole Green's function, G/ l; is con- 
trolled by a new non-zero boundary scaling dimension, 



rjh . Consequently, the hole quasiparticle residue vanishes 
at the zero temperature (T) quantum critical point due 
to an "orthogonality catastrophe" induced by a cloud of 
low energy spin excitations around the hole. The motion 
of the hole charge does not play a dominant role: this 
may be viewed as a form a "spin-charge separation" , but 
is not 2d electron fractionalization in the sense of HQ. 

As has been discussed elsewhere 0j, spin fluctuations 
near the magnetic quantum critical point in the super- 
conductors or insulators of interest are described by the 
following continuum Hamiltonian 

Utt, = J d 2 x [(tt 2 + c\(d Xl ct> a ) 2 + cl(d X2 4> a f) /2 

+ (r/2)0 2 +( 3 /24)(^ 2 ) 2 ] (1) 

Here 4> a (x, t) (a — 1, 2, 3) is the component of the magne- 
tization at wavevector G (usually G = (tt,tt)) at spatial 
co-ordinate x and time t, 7r Q is its canonically conjugate 
momentum, the only non-zero equal-time commutation 
relation is [4> a (x, t), 7173(2;', t)] = i5 a/ 35 2 (x — x 1 ), are 
velocities (possibly unequal to account for a quasi Id spa- 
tial anisotropy), r is the parameter which tunes the sys- 
tem across the quantum critical point ((4>a) 7^ in the 
magnetically ordered phase for r < r c , and (4> a ) = oth- 
erwise) , and g is the crucial quartic self interaction which 
is relevant below 3d and is responsible for the non-trivial 
universal scaling properties of the quantum critical point. 
There are no explicit damping terms for <f> a fluctuations 
in (Q) because the systems of interest do not have a gap- 
less particle- hole continuum near the wavevector G [|l0| . 

We now inject a spin S hole into the system described 
by H,/,. We consider the hole spectrum near points of 
higher symmetry in the Brillouin zone, with momentum 
k = Ki, where its dispersion has a vanishing k derivative 
i.e. near van Hove points like Ki = (it, 0). We denote the 
hole creation operator in the vicinity of Ki by ip\ a (x, t), 
with a = — S, . . . , S a spin index (for the case where the 
ground states are superconductors, yb\ a is the creation 
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operator for the Bogoliubov quasiparticles) . Coupling to 
<p a fluctuations will scatter the hole to momenta near 
K2 = Ki + G: we denote the hole creation operator near 
K2 by ipl a (x, t), and it is assumed that K2 is also a van 
Hove point. The hole Hamiltonian is then 



H h = / d 2 x 



= 1,2 V m=l,2 / 



(2) 



Here L a arc the familiar (25 + 1) x (25 + 1) angular 
momentum matrices, the a m i determine the band curva- 
tures near K1.2, and 7 is the coupling to the <p a fluctu- 
ations. Higher order couplings between the <p a and %px,2 
are also possible, but are easily shown to be irrelevant un- 
der the renormalization group (RG) discussed below. We 
assume that the hole energies, e$ satisfy t\ = £2 = £o > 
(equality holds with Ki = (tt, 0), K2 = (0, 7r), and square 
symmetry), for otherwise the hole scattering by <p a fluc- 
tuations is non-singular, and direct perturbative compu- 
tations of the hole spectrum are adequate. The singular- 
ity in the value of eo at the magnetic transition at r = r c 
is weak and subdominant, and can be safely neglected. 

A key observation follows from a simple, tree-level, RG 
analysis of Ti$ + TLh- We know that Tt^ at r = = r c is in- 
variant under the rescaling transformation x — > xe~ e , 
t — » te~ zl with 0=1. Applying this to Tth we see im- 
mediately that the a m i flow as da m i/d£ = —a m f. so the 
band curvatures are irrelevant for the low-energy theory, 
and the hole may be viewed as dispersionless. 

Before embarking on a complete RG analysis of Ti^ + 
TLh, we make some qualitative observations on the hole 
spectrum for r < r c , r > r c , and r = r c . 

For r < r c , hole motion in the magnetically ordered 
state has been studied earlier pl]| . Because (<p a ) ^ 0, 
there is a non-zero mean matrix element between the 
ipi a and ip2a states, and we have to rediagonalize TLh to 
obtain the bare hole dispersion. Simple considerations of 
energy and momentum conservation show that there is 
an infinitely sharp quasiparticle pole in the vicinity of the 
absolute band minimum: the slow quadratic dispersion 
of the hole prevents decay by emission of linearly dis- 
persing spin waves. At frequencies, to, above this pole, 
there is an incoherent gapless continuum, but its spectral 
weight vanishes rapidly as lo approaches the quasiparti- 
cle pole: after the rediagonalization of TLh required by a 
nonzero ((p a ), it is easy to see that the matrix element 
for emission of small momentum, Goldstone, spin wave 
modes is suppressed by powers of the momentum. Sim- 
ilar considerations also apply at other van Hove points 
which are not global minima. However, depending upon 
non-universal details of the band structure, in some cases 
it may be possible for the hole to emit large wavevector, 



high energy <p a quanta, and this would broaden the quasi- 
particle pole; for TLh such processes occur if one of the 
a m i < 0, and require momenta of order ci^/ct m i or larger 
(assuming TLh still applies at such momenta). Given the 
irrelevance of the a m i, the remainder of this paper will 
neglect this non-universal decay. We will only consider 
low energy <p a quanta, and assume that the high energy 
processes are either not present, or contribute a small, 
background, quasiparticle decay rate. If the latter were 
not true, there would be no sharp quasiparticle-like peak 
or threshhold in the hole spectrum, and the analysis of 
this paper would not be necessary. 

Closely related considerations apply in the spin-gap 
phase with r > r c [[l2 13 1, but with some important dif- 
ferences: (i) the incoherent continuum is separated from 
the quasiparticle pole by at least the spin gap energy, 
A ~ (r — r c ) v , where v is the correlation length exponent 
of the phase transition in TL^,; (ii) the <p a quanta are no 
longer Goldstone modes, and so the matrix element for 
emission of a (p a quantum by the hole does not vanish at 
zero momentum transfer. 

The main purpose of this paper is to understand the 
nature of the hole spectrum at the T = quantum critical 
point at r = r c , and its associated T > quantum-critical 
region. The <p a are now gapless critical excitations, but 
not Goldstone modes. Consequently, there is no factor 
of a small momentum suppressing their emission by the 
hole, and perturbative corrections in 7 are infrared sin- 
gular, as has also been noted by Sushkov JI3]. Our RG 
analysis will identify the scale-invariant quantum field 
theory which permits a resummation of the perturbative 
expansion, and shows that there is no quasiparticle pole 
at T — - and r = r c \ instead 



G h {io) = -A(e - lo) 



-1+1 Hi 



(3) 



where A is a non-universal amplitude, and the universal 
exponent rjh is computed below. 

For the critical theory, we can set a m i — 0. Further, 
examination of the perturbation theory in 7 shows that 
the presence of two hole flavors, i = 1,2, makes no mate- 
rial difference to the critical singularities: the <p a ~ -» — (p a 
symmetry of Ti^ ensures that the hole self-energy has 
terms only in even powers of 7 for which the hole flavor 
returns to its original value. Consequently we can drop 
the i index, and refer to a generic dispersionless hole ip a . 
After injection into the antiferromagnet at x = (say), 
the ip a charge will remain localized at x = 0, and its 
spin will couple to the <p a fluctuations. So we are led to 

TLa 



consider the Bose-Kondo-like model |l , Ti$ + Tts 

of a single quantum spin, S a , coupled to the bosonic <p a 
fluctuations where 



Us = "fS a (p a (x = 0), 



(4) 



[S a ,Sp] = ie a p 1 S 1 , and S a S a = 5(5+ 1). The charge 
density of the injected hole can couple only to the spin- 
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rotation invariant <pzAx — 0), and such a term is irrele- 
vant under the RG ||. 

The properties of TL^ + TLs have already been stud- 
ied in detail |J in the different physical context of Zn/Li 
impurities in the cuprate superconductors; in this ear- 
lier case the S a spin was permanently confined near the 
impurity, while in the present situation there is no impu- 
rity and S a is the spin of the injected hole. Indeed, the 
relationship between TL^ + TLs and TL§ + TLh is similar 
to that between the familiar fermionic Kondo and X-ray 
edge problems. However, the analogy is not perfect: it 
is conventional in the X-ray edge problem to neglect the 
spin of the injected hole, and merely couple its charge 
density to the fermionic bath. Here, the spin exchange 
with the (j> a quanta is paramount in both cases. 

The critical point of TL^ defines a 2+1 dimensional, 
conformally invariant field theory, and TLs is a "bound- 
ary" perturbation along the line in spacetime at x = 0. 
This perturbation flows to a fixed point which is invari- 
ant under conformal transformations which leave x = 
fixed. Correlations of S a are characterized by its bound- 
ary anomalous dimension rj'/2, which was computed to 
two- loop order in an expansion in e = 3 — d in || . 

We are interested here in the Gh = (ip a ipt), and by (||) 
we identify r/h/2 as the boundary anomalous dimension 
of ip a , which is not simply related to rj' . Gh involves an 
overlap between eigenstates of H4, + TLs and states in 
which the hole has been removed (the latter are outside 
the Hilbert space of TL^ + TLs)- The needed results do 
not follow from the previous analysis of TL$ + TLs alone, 
and require extensions we describe here. 

At one-loop order, we perform a standard momentum- 
shell RG of TLcf, + TLh in d spatial dimensions, in which 
fields with momenta between A and Ae~ e , and all fre- 
quencies, are integrated out. This is followed by the 
rescalings x — * xe~ e , t — > te~ e , <f> a — * Q e' d_1+ ' , ^/ 2 </) Q ,, 
ip a -> ^ a e {d+rnh)l/2 . The RG flow of the bulk cou- 
plings r, g is well known: g approaches a finite fixed 
point value g* at r = r c , while rj = at one 
loop order. For the hole, evaluation of the one-loop 
graphs in Fig |l| shows that rjh = S(S + 1)7 2 , where 
7 = 7(c 1 c 2 )- d/4 A- e / 2 (r(d/2)(47r)' 1 / 2 )- 1 / 2 obeys the 
same flow equation as that obtained earlier for TL$ + TLs 
[HQ (as expected): 

dj/d£ = ej/2-j 3 . (5) 

We observe that 7 also approaches a fixed-point value, 
and at which we have the anomalous dimension 

r ]h = eS(S+l)/2 + 0(e 2 ). (6) 

At the same order, a closely related computation shows 
that 77' = e ||. The extension of these results on the 
boundary exponents to two (and higher) loops requires 
the field-theoretic renormalization group and the results 



are presented elsewhere |16fl : at next order there is in- 
terference between bulk and boundary interactions, and 
the flow equation for 7 involves g. This interference is 
a novel feature of the present problem and is absent in 
the fermionic Kondo and X-ray edge problems, where the 
bulk degrees of freedom can be represented by free fields. 
Direct numerical evaluation of the two-loop corrections 
for the physical values e = 1 and 5 = 1/2 shows signifi- 
cant changes from the one loop values: % changes from 
0.375 to 0.087, while 77' changes from 1 to 0.232. 

We also computed r\h from a continuous time world- 
line quantum Monte Carlo |l7j simulation of a double 
layer Heisenberg antiferromagnet at its quantum critical 
point |l8| . We measured Gh (t) (t is imaginary time) for 
a single hole on a fixed site i, by relating it to correlators 
of spin world lines at site i which do not flip in time r. We 
obtained the estimate rjh = 0.087±0.040 from simulations 
on a 64 x 64 system by fitting to Gh (r) ~ r~ nh cxp(— e r), 
in the range 2 < Jt < 6 (J is the intralayer exchange, 
and to/ J = 2.102±0.016 was determined separately from 
the ground state energies of the antiferromagnet with and 
without a hole) ; the numerical data is in Fig ||. 

In the r > r c , T > vicinity of the critical point, (||) 
generalizes to 

GhM - 5^** (^A (7) 

where A and en are the same constants appearing in 
(they have absorbed a non-universal renormalization 
from the coupling to (ft a modes) , while <&h is a completely 
universal function. We obtained numerical results for $^ 
using the large N method discussed in |)| and the results 
are shown in Fig ||. The k dependence of Gh arises only 
from the irrelevant a m i couplings, and their main effect 
is to replace en by the actual hole dispersion near the 
van-Hove point. 

We conclude by discussing possible physical applica- 
tion of our results to photoemission measurements on 
the cuprate superconductors. We have already argued 
elsewhere ||,|o| that (P should be a reasonable descrip- 
tion of the antiferromagnetic fluctuations in the low and 
moderate doping regime, both above and below the su- 
perconducting T c . There is evidence from NMR |2l]] and 
neutron scattering experiments [ p2[ that the quantum 
critical region of the r = r c critical point describes the 
antiferromagnetic fluctuations above T c ; so for the same 
systems, the r = r c spectrum in Fig || should apply to 
photoemission at the anti-nodal points. The proximity 
of such a magnetic quantum critical point may be asso- 
ciated with the onset of quasi Id correlations, but this 
incidental to our theory — the anisotropy is merely re- 
flected by changes to the couplings in TL<$,. Below T c , 
the measured anti-nodal spectrum p3| is similar to the 
r > r c spectrum in Fig |^: this is accounted for in our 
approach by the reasonable assumption that the onset of 
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superconductivity induces the spin-gap-like correlations 
and so increases the value of the effective r controlling 
the magnetic fluctuations. The high frequency tail of the 
EDC both above and below T c should decay as l/ui 1 ~ Vh . 
Note that this tail is present for any non-zero value of r/h, 
with a non-universal amplitude determined by A. How- 
ever, our estimate of the value of rjh here is quite small, 
and it appears that the interaction effects discussed here 
cannot explain the prominent tail seen in current exper- 
iments on their own. It is likely that the strong disorder 
in the local gap values (as seen in recent STM measure- 
ments p4|) also plays an important role. 
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was supported by NSF Grant DMR 96-23181 (USA), the 
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FIG. 1. One- loop Feynman diagrams. The full line is tp a 
and the dashed line is d> a . 
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FIG. 2. Quantum Monte Carlo results for a localized hole 
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layer Heisenberg antiferromagnet. Both axes are logarithmic. 
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FIG. 3. Large N Q results for the imaginary part of 
the scaling function in (0). At large N rjh — 1/2. 
We used an interpolation form for the <f> a propagator 
l/(k 2 — uj 2 + m 2 — 2iFuj), where m/T and T /T are univer- 
sal functions of A/T Both m and the damping, V, are 
non-zero even at r — r c because the thermally excited <j> a 
quanta scatter strongly off each other by the non-zero fixed 
point value of the bulk interaction, g, in Ji^. For r = r c we 
used m = T = T, and for r > r c we used F = m/5 = T. 



5 



